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We present a study of the local strain effects associated with vacancy defects in strontium titanate 
and report the first calculations of elastic dipole tensors and chemical strains for point defects in 
perovskites. The combination of local and long-range results will enable determination of x-ray 
scattering signatures that can be compared with experiments. We find that the oxygen vacancy 
possesses a special property — a highly anisotropic elastic dipole tensor which almost vanishes 
upon averaging over all possible defect orientations. Moreover, through direct comparison with 
experimental measurements of chemical strain, we place constraints on the possible defects present in 
oxygen-poor strontium titanate and introduce a conjecture regarding the nature of the predominant 
defect in strontium-poor stoichiometrics in samples grown via pulsed laser deposition. Finally, 
during the review process, we learned of recent experimental data, from strontium titanate films 
deposited via molecular-beam epitaxy, that show good agreement with our calculated value of the 
chemical strain associated with strontium vacancies. 

PACS numbers: 61.72.Bb,61.72.Hh,61.72.jd,62.20.D- 



I. INTRODUCTION 

Perovskites in general, and strontium titanate in par- 
ticular, are some of the most frequently and exhaus- 
tively studied materials in solid-state physics and chem- 
istry This attention has largely derived from their di- 
verse and interesting pr ope rties: high piezoelectricity,^ 
quantum paraelectricitypEl ferroelectr icity ffl^ uniaxial 
stress,^ and colossal magnetoresistanceP^ Further, the 
cubic perovskites manifest intriguing effects of underly- 
ing quantum fluctuations, since, although they are com- 
prised of relatively heavy atomic constituents, a number 
of competing structures are energetically and structurally 
similar.^] Characterization of the low-temperature order 
parameters of these materials remains an open and en- 
gaging question! 1 ^ 

Strontium titanate is a model perovskite — commonly 
available and reflective of many of the above properties 
of that material family. Specifically while strontium ti- 
tanate is a wide band-gap insulator at room temperature, 
it exhibits semiconductivity at elevated tempera tures 
through doping or non-stoichiometric compositio n 1 13 ! 14 -^ 
and superconductivity at low temp eratures through re- 
duction via addition of oxygen ESliSI The structural phase 
diagram of strontium titanate comprises a high-tempera- 
ture cubic phase and a low-temperature tetragonal phase, 
with a transition temperature near 105 K.Q2HI2] The cu- 
bic perovskite structure is particularly interesting due to 
the richness of its phase diagram (nonpolar antiferrodis- 
tortive to ferroelectric to antiferroelectric phases) J2212U to 
the capacity of these phases to emerge from miniscule 
deviations from the cubic lattice and its skeleton of oc- 
tahedral oxygens (often through rigid rotations of such), 
and to open questions regarding the types (displacive or 
order-disorder) of the transitions among these phases 



Defects and vacancies play a particularly important 
role in the chemistry of perovskites and deserve contin- 
ued study in strontium titanate due to the electronic and 
superconducting effects of doping as well as their role in 
the interface region of heterostructures.^l In the dilute 
limit, the mechanics of defects are fully determined by an 
examination of stress-strain effects, in particular the elas- 
tic dipole tensor, which motivates this work's emphasis 
on such a quantity. Our presentation of both short-range 
displacements around a point defect as well as long-range 
effects (characterized by the elastic dipole tensor) allows 
for the calculation of x-ray scattering signatures. These 
quantities also enable the prediction of defect mechan- 
ics, such as the behavior of defects within an externally 
imposed strain gradient (as present in heterostructures) , 
as well as the ratio of chemical strain to stoichiomet- 
ric deviation (a direct experimental observable) . Finally, 
through comparison of our predictions of chemical strain 
with experimental results, we draw a number of impor- 
tant conclusions regarding the nature of point defects in 
non-stoichiometric strontium titanate. 



II. BACKGROUND 

Point defects introduce lattice distortions on both local 
and long-range scales. While the short-range distortions 
must be described by a potentially large set of atomic 
displacements, the long-range elastic distortions may be 
completely described by a single tensor, the elastic dipole 
tensor.^ 

The clastic dipole tensor and its relation to clastic ef- 
fects may be understood by the following simple consid- 
erations. Consider to quadratic order the most general 
expansion of the free energy per unit volume of a crys- 
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tal in terms of the strain (i and j refer to coordinate 
axes) and the number of defects per unit volume rid, 
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The Taylor expansion coefficients Ed, Edd, C%jki, an d 
Gij are, respectively, the defect formation energy, an av- 
erage inter-defect interaction energy, the components of 
the elastic stiffness tensor of the material, and the com- 
ponents of the elastic dipole tensor of the defects. The 
negative derivative of the free energy with respect to 
strain then gives the stress, 

— o — = a ij — ~ ^ Cijki tki + nd Gij. (2) 

To isolate the elastic dipole tensor, we can consider the 
rate of change of the stress in the crystal per unit concen- 
tration of defects, while holding strain fixed, that is, un- 
der strain control. Although challenging experimentally, 
strain control is quite convenient computationally since 
it corresponds to performing calculations with fixed lat- 
tice vectors. This derivative thus gives the elastic dipole 
tensor G directly, 

da. 



dn d 
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so that positive diagonal components of G indicate that 
the presence of defects tends to expand the crystal along 
the corresponding directions. 

Alternatively, we can also consider the derivative of 
the strain in the crystal per unit concentration of defects 
under stress or load control (holding stress fixed). While 
stress control is computationally more complicated than 
strain control, it is the most common experimental situa- 
tion. Under experimental conditions, the crystalline lat- 
tice vectors relax such that there is essentially zero stress 
(under normal laboratory conditions, atmospheric pres- 
sure corresponds to a negligible stress). This criterion 
allows the relation of strain to a newly defined quantity, 
A, 



dei 



dn d 



^ Sijki Gki = Ai 



(4) 



ki. 



where S, the elastic compliance tensor, is the inverse of 
the elastic stiffness tensor C, and A is defined as the 
defect-strain tensor, which is the strain per unit defect 
concentration induced in a crystal at fixed stress. 

Relying upon the above derivations, the numerical cal- 
culation of the clastic dipole tensor G is straightforward. 
We compute the stress induced with the introduction of a 
single defect in a supercell (maintaining fixed lattice vec- 
tors, but allowing relaxation of the atomic coordinates). 
From Equation pj, this yields 

*«) (5) 

V A o~a , 



where a d and a are the stresses with and without the 
defect in the cell, respectively, and V is the supercell 
volume. Once G is known, A may also be computed 
directly from Equation Q. As a practical matter, we 
note that in this approach, the lattice vectors need not be 
those of a fully relaxed bulk crystal, provided the strain 
is small and kept fixed. 

Experimental works often report the variations in 
chemical strain, the strain due to the presence of de- 
fects, with respect to stoichiometric deviations in the 
crystalline chemical formula. From the above consider- 
ations, the chemical strain is e c = n^A. Deviations in 
stoichiometry specify the number of defects per chemi- 
cal unit, 5, so that, in this context, the concentration of 
defects per unit volume is nd = S/V c , where V c is the 
volume of the chemical unit. These two relations then 
immediately provide the chemical strain as proportional 
to this stoichiometric defect deviation, 5, 
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Experimentally, one does not generally obtain a full ten- 
sor for e c but, instead, an average over all equivalent 
defect orientations which restore the symmetry of the un- 
derlying crystal. In cubic crystals, one measures a scalar 
e c which corresponds to the mean diagonal component of 



III. METHODS 

To simulate strontium titanate, we employ a shell-po- 
tential modeP^ parameterized for strontium titanate.^ 
Shell-potential models are formulated as an extension to 
ionic pair potentials and employed to capture the po- 
larizability of the atomic constituents. The shell model 
separates each ion into two parts, a core and an outer 
shell, which possess individual charges that sum to the 
nominal charge of the ion. The total model potential U 
consists of three terms, 



U = Up + U c 



(7) 



representing, respectively, the polarizability of the ions 
and the Coulomb and short-range interactions among the 
ions. The polarizability is captured by harmonic springs 
connecting the core and shell of each ion, so that Up has 
the form, 



Up 
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where |Arj| is the core-shell separation for ion i and the 
ki are a set of ion-specific spring constants. The Coulomb 
contributions take the form, 
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Shell Core Spring Constant 

Charge [e] Charge [e] [eV • A~ 2 ] 



Sr 2+ 1.526 0.474 
Ti 4+ -35.863 39.863 
2 ~ -2.389 0.389 


11.406 
65974.0 
18.41 


TABLE I: Electrostatic shell-model potential parameters for 
strontium titanate (from Akhtar et al.^). 


Interaction A [eV] p [A] 


C [eV • A 6 ] 


Sr 2 + ^ O 2 - 776.84 0.35867 
Ti 4+ ^ O 2 - 877.20 0.38096 
O 2 " ^ 2 ~ 22764.3 0.1490 


0.0 
9.0 
43.0 



TABLE II: Short-range shell-model potential parameters for 
strontium titanate (from Akhtar et alP^). 



where i and j range over all cores and shells (exclud- 
ing terms where i and j refer to the same ion), qi and 
qj are the corresponding charges, is the distance be- 
tween the charge centers, and k c is Coulomb's constant. 
We compute this Coulombic interaction^ using a Parti- 
cle Mesh Ewald algorith m^" 29 ! with all real-space pair- 
potential terms computed out to a fixed cutoff distance 
using neighbor tables. Finally, the short-range interac- 
tions are included through a sum of Buckingham^ pair 
potentials (which can be viewed as combinations of Born- 
Mayer^ and Lennard-Jones^l potentials) of the form, 

u b = \ E (Aije-^ - ( ;,,,;;) , (10) 
id 

where i and j range over all shells and A^j, pij, and 
are pair-specific adjustable parameters. Here, the first 
term (Born-Mayer) serves as a repulsive short-range in- 
teraction to respect the Pauli exclusion principle, and 
the second term (Lennard-Jones) models the dispersion 
or van der Waals interactions.^ The specific electro- 
static and short-range shell-model parameters used in 
this study were fit to strontium titanate by Akhtar et 
al.pSl with values as listed in Tables [i] and [n] Finally, 
we wish to emphasize again, as it is rarely mentioned 
explicitly in the shell-potential literature, that the pair- 
potential terms in Ub apply to the shells only, and not 
to the cores. 

Shell models have been extensively used for decades as 
the primary empirical potential for modeling perovskites 
and other oxides JHHU We tested the correctness of our 
coded implementation of this potential through compar- 
isons of lattice constants and elastic moduli and find ex- 
cellent agreement. For instance, using the same shell 
potential and ground-state structure, we predict a lat- 
tice constant for cubic strontium titanate of 3.881 A, 
which is within 0.3% of the value calculated by Akhtar 
et alPFor the elastic moduli, we find C u = 306.9 GPa, 
C12 = 138.7 GPa, and C 44 = 138.8 GPa, which are 



within 1.8%, 1.0%, and 0.7%, respectively, of the val- 
ues from Akhtar et alP^From this, we conclude that our 
implementation of the potential is correct. 

We further note that the static dielectric constant of 
216.99, as calculated by Akhtar et al.^for the same pa- 
rameter set as our shell potential, is 28% lower than the 
experimental value for strontium titanate of 301.00, a 
level of agreement typic al of results from empirical po- 
tentials. As GillarpSGU and StonehanPH discuss, there 
exists an important electrostriction effect whereby the 
tendency of the crystal to screen electric fields impacts 
the elastic dipole tensor. (Note that this effect scales 
as 1 — e _1 , where e is the static dielectric constantP^) 
Given the high dielectric constant of strontium titanate, 
this screening is nearly perfect in both the actual exper- 
iment and our model case; therefore, we expect that this 
effect is captured well in our calculations below, despite 
our relatively large fractional error in the static dielectric 
constant. 

As is well known, strontium titanate has a large num- 
ber of similar, competing ground-state structures.^ We 
should emphasize, at this point, that the main quantities 
of interest to this study, either local atomic displacements 
or the clastic dipole and defect-strain tensors [from Equa- 
tions |3]H([5])], are all defined as defect-induced changes 
relative to the bulk structure and so are likely quite in- 
sensitive to which of the competing structures are used 
to represent the bulk. 

Accordingly, we have carried out what we regard as 
a thorough, but not exhaustive, search for a probable 
ground-state structure. Indeed, we have found no al- 
ternative structure which relaxes to an energy less than 
our candidate ground-state structure within our poten- 
tial. We performed quenches on hundreds of random dis- 
placements from the idealized positions of the lxlxl 
primitive unit cell to explore various potential reconstruc- 
tions for supercells up to 6 x 6 x 6. We also consid- 
ered a number of highly ordered, human-inspired config- 
urations commensurate with the antiferr odisto rtive dis- 
ordering that is o bserved experimentally*^ 9 -* 4 ^ and pre- 
dicted theoretically! 20 * 41 -^ Among those minima which 
we explored, we selected the lowest-energy configura- 
tion to serve as the bulk crystalline state throughout 
this study. This configuration possesses a fairly regular 
pattern, namely, each oxygen octahedron rotates slightly 
along ±[111] (trigonal) directions in an alternating three- 
dimensional 2x2x2 checkerboard pattern. To aid in vi- 
sualizing this reconstruction, Figure [l] depicts the atomic 
structure of reference non-reconstructed cubic strontium 
titanate, with the oxygen octahedra lattice indicative of 
perovskite materials. 

We investigated five defects in strontium titanate: the 
oxygen, strontium, and titanium vacancies and the stron- 
tium-oxygen and titanium-oxygen divacancies. Since the 
octahedra rotations break the original crystalline sym- 
metry and generate a set of different symmetry-related 
reconstructions, each of these five defects can be situated 
in multiple equivalent sites within each reconstruction. 
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FIGURE 1 


Atomic 


structure of ] 


ion-reconstructed cubic 



strontium titanate: oxygen octahedra (^) surrounding ti- 
tanium atoms (♦), with strontium atoms (■) outside the 
octahedra and equidistant from the titanium sites. 



Even well below room temperature, strontium titanate 
shows local fluctuations among the possible reconstruc- 
tions. Thus, in addition to raw results for a specific real- 
ization of each defect in a given reconstruction, we also 
report results for each defect averaged over all possible 
reconstructions for a given orientation of the defect. The 
strontium and titanium vacancies do not select a spe- 
cific direction and, thus, this averaging represents the full 
cubic crystalline symmetry group; their respective ten- 
sors therefore are always diagonal with cubic symmetry. 
The other defects do select specific crystalline directions, 
which must be specified when reporting the correspond- 
ing defect tensors. We thus reconstructionally average 
these latter defect tensors using either symmetry argu- 
ments or explicit calculations with different reconstruc- 
tions, as appropriate. Finally, with a view to chemical 
strain measurements in macroscopic samples, we also re- 
port final averages over all defect orientations. 

In the case of the oxygen vacancy, the oxygen sits be- 
tween two unique nearest-neighbor titanium atoms, thus 
uniquely distinguishing the Ti-Vo-Ti direction, which we 
define as [100], among the three cubic axes. Next, the 
titanium-oxygen divacancy selects a unique Vti-Vo direc- 
tion, which we define as [100], directed from the titanium 
site toward the oxygen site. Finally, the strontium-oxy- 
gen divacancy selects a unique Vs r -Vo direction, which 
we define as [110] , directed from the strontium toward the 
oxygen site. Once the reconstructional averaging is ac- 
complished, the average over defect orientations requires 
generation of the crystal's response to all different pos- 
sible orientations (three for the oxygen vacancies, six for 



the titanium-oxygen divacancy, and twelve for the stron- 
tium-oxygen divacancy) and restores full cubic symme- 
try. 

To obtain the results in Section |IV| each of these 
defects was placed within the bulk-reconstructed stron- 
tium titanate supercell, with cubic symmetry as exper- 
imentally observed for T > 105 K,E2Hl9] anc j then re- 
laxed via the technique of preconditioned conjugate gra- 
dient minimization^^] (specifically, the Polak-Ribiere^ 
method) to find the minimal energy configuration (to 
within ~1 |^eV). Supercell convergence studies examined 
all five defects in cells containing up to 13 720 atoms 
and verified that such defects were sufficiently separated 
so that interactions across periodic supercell boundaries 
were negligible.^ The final relaxed atomic configurations 
in the largest cells (14 x 14 x 14) provide the local strain 
fields which we report below for each defect. To deter- 
mine the long-range strain fields, we compute the elas- 
tic dipole tensor G through the prescription described 
in Equation ^ above; namely, we calculate the stress 
induced by the introduction of a single defect in the su- 
percell, holding the lattice vectors fixed while allowing 
the atomic coordinates to relax. 



IV. RESULTS 

For a series of important and fundamental strontium 
titanate defects, we examine both the elastic dipole and 
defect-strain tensors, as introduced in Section [TTJ as well 
as the local strains surrounding each defect. This section 
first examines the role of the oxygen vacancy as a case 
study of a defect in strontium titanate. Subsequently, 
the results for the same set of studies are presented for 
four other point defects in strontium titanate: strontium 
and titanium vacancies and strontium-oxygen and tita- 
nium-oxygen divacancies. Section [V] continues with more 
general implications of our results. 

A. Oxygen vacancy 

As described above, there are three distinct orienta- 
tions for the oxygen vacancy as defined by the Ti-Vo-Ti 
direction. Moreover, because of the reconstruction, there 
are in fact two distinct classes of site within each possible 
orientation. As one moves in the positive sense of direc- 
tion along the Ti-Vo-Ti axis, these sites are distinguished 
by whether the rotation of the octahedra surrounding the 
two titanium sites changes from positive to negative or 
from negative to positive (in the right-hand sense about 
the +[111] direction). Below, we report results for the 
latter type of site. 

We first examine the elastic dipole tensor computed 
according to Equation To explore convergence, we 
compute the dipole tensor components in supercells of 
sizes 2x2x2, 4x4x4, ... , 14 x 14 x 14, contain- 
ing between 40 and 13 720 atoms with defect separations 
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FIGURE 2: Diagonal components of elastic dipole tensor for 
oxygen vacancy in strontium titanate. Data show linear con- 
vergence with inverse linear dimension of supercell size. 



between ^8 and ~54 A. 

Figure [2] depicts the convergence of the diagonal com- 
ponents of the elastic dipole tensor as a function of in- 
verse linear dimension of the supercell. The linear behav- 
ior in the figure for large cells indicates that this quan- 
tity converges in the same way as the Coulomb inter- 
action between defects and allows extrapolation of the 
fully converged value for these components in an infinite 
supercell. We observe exactly the same linear behavior 
with inverse linear dimension of cell for the convergence 
of the off-diagonal components of the dipole tensor. A 
linear fit to the data, for all tensor components for cell 
sizes in the range exhibiting linear behavior (6x6x6 
through 14 x 14 x 14), yields the elastic dipole tensor, 
extrapolated to infinite cell size, for this oxygen vacancy 
in strontium titanate, 



4.53 -3.11 -3.11 
-3.11 -2.13 1.06 1 < Y. 
-3.11 1.06 -2.13 



Go = 



As described above in Section at finite tempera- 
tures (T > 105 K), local fluctuations in the reconstruc- 
tion make it appropriate to average this tensor over all 
reconstructions. The result is to eliminate the off-diago- 
nal components, leaving the reconstructionally averaged 
elastic dipole tensor, 
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Since the underlying, non-defected, crystal is now cubic, 
we can readily apply Equation Q to Go to obtain the 
reconstructionally averaged defect-strain tensor, 



16.33 0.00 0.00 
0.00 -8.05 0.00 
0.00 0.00 -8.05 



• 




• 




• 




• 


♦ 


\J 


♦ 


o 


♦ 


U 


♦ 


{J 




















r\ 
vJ 




r\ 
w 












ft* 








o 




O 




o 




• 




• 




• 




• 



FIGURE 3: Local strain pattern for oxygen vacancy in the 
eii.2 plane of strontium titanate: titanium atoms (♦), oxy- 
gen atoms (#), oxygen vacancy ( )). Atomic displacements 
exaggerated by a factor of three for clarity (-»■), displayed for 
significant in-plane displacements (> 0.1 A) only. 



The above results indicate that the oxygen vacancy 
tends to cause the crystal to expand along the Ti-Vo-Ti 
direction and to contract along the two orthogonal di- 
rections by an amount which results in negligible net 
volume change in the crystal. When the above tensor 
is averaged over all defect orientations (permutations of 
the three coordinate axes) , the result is near perfect can- 
cellation, resulting in a constant tensor (multiple of the 
identity) with a uniform chemical strain per unit defect 
concentration of +0.07 A 3 . This corresponds to a ratio 
of chemical strain e c to the deviation from stoichiometry 
S in SrTi03_a of e c /S — +0.001, indicating a very slight 
tendency for the crystal to expand due to the presence 
of oxygen vacancies. 

Following Figure [3j we now examine the local strain 
of this system after reconstructional averaging. We note 
that the removal of the oxygen ion, with a nominal charge 
of —2, leaves an effective local positive environment in the 
location of the vacancy. We should expect pronounced 
Coulombic response to this environment in the form of lo- 
cal crystal polarization. Indeed, the nearest neighbors of 
the oxygen vacancy, the two titanium atoms, now move 
directly away from the vacancy on the precise vector con- 
necting them, by 0.21 A. The next-nearest neighbors are 
the eight oxygens comprised of two equidistant rings of 
four oxygens each, each of which move 0.21 A toward 
the vacancy (0.19 A along e\ with remaining projection 
of 0.08 A either along e 2 or e 3 as dictated by symme- 
try), and one ring of four strontium atoms, each of which 
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FIGURE 4: Diagonal components of elastic dipole tensor for 
strontium vacancy in strontium titanate. Data show linear 
convergence with inverse linear dimension of supercell size. 
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moves 0.10 A directly away from the vacancy. The third 
shell of neighbors, a set of six oxygens, each one lat- 
tice constant away from the vacancy along the three lat- 
tice directions in the crystal, moves different amounts 
depending upon the vector; the two oxygens along the 
Ti-Vo-Ti direction move 0.06 A directly away from the 
vacancy (pushed sterically by the second-shell titanium 
atoms), while the other four oxygens move by a negli- 
gible amount (only 0.01 A). While the fourth shell of 
neighbors, a set of eight titanium atoms at a distance 
v5/2 a from the vacancy, moves negligibly, the fifth shell 
of neighbors (sixteen oxygen and eight equidistant stron- 
tium atoms) shows significant movement in the strontium 
atoms of 0.15 A away from the vacancy (0.08 A along 
ex with remaining projections of 0.09 A along both e% 
and £3), even though the oxygen atoms move a negligi- 
ble amount. Finally, in the sixth shell of neighbors, a 
set of twelve oxygen atoms, eight oxygens (those not in 
the plane containing the defect and perpendicular to e\) 
move 0.10 A toward the vacancy (0.02 A along e\ with 
remaining projection of 0.10 A either along e-i or £3 as 
dictated by the symmetry), while four others (those in 
the plane perpendicular to e{) move a negligible amount. 
All other atoms in the crystal move less than 0.06 A. 



Finally, we would like to comment on the relation be- 
tween local strains and defect tensors. We observe that 
the direction of motion of the near-neighboring atoms of- 
ten correlates with the far-field motion described by the 
defect-strain tensor. In this case, in the first and second 
shells, we see a general pattern of movement which is 
away from the defect along e.\ and toward the defect in 
the other two directions, consistent with the signs in the 
long-range defect-strain tensor. 



FIGURE 5: Local strain pattern for strontium vacancy in 
the eie2 plane of strontium titanate: strontium atoms (□), 
oxygen atoms (•), strontium vacancy (□). Atomic displace- 
ments exaggerated by a factor of three for clarity (->■), dis- 
played for significant in-plane displacements (> 0.1 A) only. 



B. Strontium vacancy 

We now repeat the above procedures to obtain similar 
results for the strontium vacancy. As described above, 
the strontium-vacancy site defines no unique direction 
and reconstructional averaging recovers the full cubic 
symmetry group. For any realization of the reconstruc- 
tion there are in fact two distinct types of strontium sites. 
Each such site sits at the center of a cube with oxygen oc- 
tahedra at the corners with alternating signs of rotations. 
The results reported below, prior to reconstructional av- 
eraging, correspond to the site in which the rotation at 
the [111] corner is positive. 

For the elastic dipole tensor we find 
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where Figure [4] shows the convergence of the diagonal 
elements of the above tensor. Performing the reconstruc- 
tional average gives 



G Sr = 
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with a corresponding defect-strain tensor, 
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FIGURE 6: Diagonal components of elastic dipole tensor for 
titanium vacancy in strontium titanate. Data show linear 
convergence with inverse linear dimension of supercell size. 



This result expresses the tendency of the crystal to ex- 
pand due to the strontium vacancy, by an amount signif- 
icantly greater than the net effect of the oxygen vacancy. 
Since the original defect defines no unique direction, no 
orientational averaging is necessary, and we find a ratio 
of chemical strain e c to the deviation from stoichiometry 
5 in Sri-aTiOa of e c /8 = +0.030, indicating a tendency 
for the crystal to expand due to the presence of strontium 
vacancies. 

We now examine the local strains around such a stron- 
tium vacancy after reconstructional averaging (Figure pp. 
Twelve oxygen atoms (three in each of four neighboring 
octahedra) are nearest neighbors to the strontium va- 
cancy; these twelve atoms all move directly away from 
the vacancy by a distance of 0.15 A. The next-nearest 
neighbors are the eight titanium atoms in each of the 
surrounding octahedra, each of which moves 0.09 A di- 
rectly toward the strontium vacancy. All other atoms 
move less than 0.08 A. 

In this expected, both the defect-strain ten- 

sor and the local strain displacements are symmetric af- 
ter reconstructional averaging. The nearest-neighbor dis- 
placement shows an expansion in all directions, similar 
to the defect-strain tensor. The next-nearest neighbors, 
which have an opposite charge from that of the near- 
est neighbors, move in the opposite direction, reinforcing 
that there is no simple connection between local displace- 
ments and far-field strain patterns. 



C. Titanium vacancy 

The titanium- vacancy site also defines no unique direc- 
tion in the ideal crystal, and the reconstructional averag- 
ing restores the full cubic symmetry. For any realization 
of the reconstruction there are, in fact, two distinct types 
of titanium site. Each such site sits at the center of an 
octahedron with either positive or negative signs of the 
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FIGURE 7: Local strain pattern for titanium vacancy in the 
ei&2 plane of strontium titanate: titanium atoms (♦), oxy- 
gen atoms (•), titanium vacancy ( ). Atomic displacements 
exaggerated by a factor of three for clarity (-*■), displayed for 
significant in-plane displacements (> 0.1 A) only. 



rotations relative to the +[111] axis. The results reported 
below, prior to reconstructional averaging, correspond to 
the site in which the rotation is positive. 

We first report our results for the elastic dipole tensor, 



Gxi = 




-0.70 
28.08 
-0.70 




where Figure [6] shows the convergence of the diagonal el- 
ements of the above tensor. The reconstructional average 
is then 



G i 



28.08 
0.00 
0.00 



with defect-strain tensor, 



A T i = 



23.92 
0.00 
0.00 



0.00 
28.08 
0.00 



0.00 
23.92 
0.00 



0.00 
0.00 
28.08 



0.00 
0.00 
23.92 



eV, 



A . 



This defect-strain tensor expresses the tendency of the 
crystal to expand due to the titanium vacancy — a sig- 
nificantly greater amount even than that of the stron- 
tium. Again, since the original defect defines no unique 
direction, no orientational averaging over defect types is 
necessary. Finally, we report a ratio of chemical strain 
e c to the deviation from stoichiometry S in SrTi^^-Oa 
of e c /5 = +0.402, indicating a significant tendency for 
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the crystal to expand due to the presence of titanium 
vacancies. 

As depicted in Figure [7J we now describe the local 
strain effects on the atoms surrounding the titanium va- 
cancy. The nearest neighbors are six surrounding oxy- 
gen atoms which move 0.22 A directly away from the 
titanium vacancy. The next-nearest neighbors are the 
eight surrounding strontium atoms along the body diag- 
onals from the titanium vacancies; these strontium atoms 
move the very significant distance of 0.52 A directly to- 
ward the titanium vacancy. The third shell is made up of 
six titanium atoms separated by a lattice constant from 
the vacancy along all three directions (positive and neg- 
ative); all of these titanium atoms move 0.09 A directly 
away from the titanium vacancy. Finally, the fourth shell 
of atoms is comprised of twenty-four oxygen atoms, ar- 
ranged as six groups of four oxygens, each in a diamond- 
shape with its center one lattice coordinate away from the 
titanium vacancy in each lattice direction. These oxygens 
each move 0.08 A away from the vacancy (0.08 A along 
the vector separating the diamond-group from the va- 
cancy, and 0.03 A along either of the two other directions, 
so as to cause the diamond- group to spread outward). All 
other atoms in the crystal move less than 0.07 A. 

We again observe connections between the reconstruc- 
tionally averaged local displacements and the far-field de- 
fect-strain tensor which have similar symmetry. In this 
case of the titanium vacancy, the nearest-neighbor atoms 
are displaced away from the vacancy, showing the same 
behavior as the defect-strain tensor. Interestingly, how- 
ever, the next-nearest neighbors, which move toward the 
titanium vacancy, actually have more than twice the dis- 
placement of the nearest neighbors. So here we observe 
that the far-field strain does not correlate with the largest 
magnitude displacement, but instead with that of the 
nearest-neighbor atoms. 



D. Strontium-oxygen divacancy 



is then 



As described above in Section pTI] there are twelve dis- 
tinct orientations for the strontium-oxygen divacancy as 
defined by the Vs r -Vo direction. For any realization of 
the reconstruction there are in fact two distinct types of 
strontium-oxygen sites. Each such strontium site sits at 
the center of a cube with oxygen octahedra at the corners 
with alternating signs of rotations. The results reported 
below, prior to reconstructional averaging, correspond to 
the [110] defect with negative sense of rotation for the 
octahedron at the [111] corner. 

The elastic dipole tensor for this defect is 



Gsro — 



-4.62 
3.00 

-2.28 



3.00 
-4.62 
-2.28 




G 



SrO 



-3.00 
1.78 
0.00 



with defect-strain tensor, 



k-SrO 




1.78 
-3.00 
0.00 



2.14 
-9.37 
0.00 





where Figure [8] shows the convergence of the diagonal el- 
ements of the above tensor. The reconstructional average 



This defect-strain tensor expresses a slight tendency of 
the crystal to contract due to the strontium-oxygen di- 
vacancy. The lower symmetry of this defect, with its ori- 
entation along a diagonal, leads to remaining off-diagonal 
elements even after reconstructional averaging. However, 
when the above tensor is averaged over all twelve defect 
orientations, the result is a constant tensor with a uni- 
form chemical strain per unit concentration of defect of 
—0.50 A 3 . This corresponds to a ratio of chemical strain 
e c to the deviation from stoichiometry 8 in Sri_,5Ti03_5 
of e c /S — —0.008, indicating a tendency for the crystal 
to contract due to the presence of strontium-oxygen di- 
vacancies. 

Referring to Figure [9] we now examine the local dis- 
placements after reconstructional averaging. The situa- 
tion with this divacancy is more complicated than that 
of earlier isolated atomic vacancies, leading us to char- 
acterize the atomic displacements with respect to each 
independent missing atom in the strontium-oxygen diva- 
cancy. The first shell is comprised of the eleven remaining 
oxygen atoms that are nearest neighbors of the strontium 
vacancy. This strontium vacancy would normally have 
twelve neighboring oxygen atoms forming three squares 
in the three planes, each comprised of four atoms cen- 
tered around the strontium vacancy; however, one of 
these oxygen atoms is missing to form the divacancy. The 
four oxygens, which are farthest from the oxygen vacancy 
but not in the same square as the oxygen vacancy, all 
move away from the strontium vacancy by 0.14 A (0.10 A 
along £3, and 0.09 A along either ei or e 2 , depending 
upon symmetry, with 0.004 A along the opposite vector, 
chosen such that most of each displacement is within the 
plane of squares to which that oxygen belongs). The four 
oxygens that are closest to the oxygen vacancy and still 
not in the same square as the oxygen vacancy also move 
as a group. These oxygen atoms all move away from 
the strontium vacancy by 0.11 A (0.05 A along (§3, and 
0.09 A along either i\ or §2, depending upon symme- 
try, with 0.05 A along the opposite vector, chosen such 
that more of each displacement is within the plane of the 
square to which that oxygen belongs). The two oxygens, 
in the same square as the oxygen vacancy and closer to 
such vacancy, are displaced away from the defect, in that 
plane, by 0.12 A (0.12 A along either e\ or e 2 , and 0.04 A 
along the opposite vector, chosen to ensure that the dis- 
placement maximizes its overall movement away from the 
remaining oxygen in this square) . The final oxygen atom 
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FIGURE 8: Diagonal components of elastic dipole tensor for 
strontium-oxygen divacancy in strontium titanate. Data show 
linear convergence with inverse linear dimension of supercell 
size. 



is also in the same square as the oxygen vacancy, but is 
the farthest atom from the vacancy in this square; it is 
displaced directly away from the defect in that plane by 
0.14 A. Also included in the first shell are the remaining 
oxygen atoms that are nearest neighbors to the missing 
oxygen vacancy. That oxygen vacancy, part of two oxy- 
gen octahedra, has eight neighboring oxygens, four in 
each octahedron; of those eight oxygens, the four closest 
to the strontium vacancy have already been considered 
as part of the nearest neighbors to the missing strontium 
atom. The remaining four oxygen atoms in the other oc- 
tahedron, farthest from the strontium vacancy, move the 
appreciable distance of 0.33 A toward the vector connect- 
ing the strontium and oxygen vacancies (0.30 A along 63, 
and 0.11 A along either e\ or e.2, depending upon sym- 
metry, with 0.005 A along the opposite vector, chosen 
such that most of each displacement in the eif§2 plane is 
toward the oxygen vacancy). Finally, this first shell also 
includes the three strontium atoms (the fourth is itself 
missing) nearest the oxygen vacancy. The one strontium 
atom farthest from the strontium vacancy moves directly 
toward the oxygen vacancy by 0.03 A (0.02 A along both 
&i and £2). The other two strontium atoms move away 
from the oxygen vacancy by 0.10 A (0.10 A along either 
ei or e2, depending upon symmetry, with 0.04 A along 
the opposite vector, chosen such that the displacement 
maximizes the distances of these strontium atoms from 
the strontium vacancy). 

The second shell is comprised of the eight titanium 
atoms that are the nearest neighbors to the missing stron- 
tium atom. Two of these nearest neighboring titanium 
atoms, which could alternatively have been categorized 
as the nearest-neighboring shell of atoms from the oxygen 
vacancy, move away from the strontium vacancy (as well 
as the oxygen vacancy) by 0.19 A (0.19 A away from the 
vacancies along £3, and 0.03 A toward the strontium va- 
cancy along both e± and e.2)- Two other nearest-neighbor 
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FIGURE 9: Local strain pattern for strontium-oxygen di- 
vacancy in the iii-2 plane of strontium titanate: strontium 
atoms (□), oxygen atoms (•), strontium vacancy (□), oxy- 
gen vacancy (O). Atomic displacements exaggerated by a 
factor of three for clarity (-»•), displayed for significant in- 
plane displacements (> 0.1 A) only. 



titanium atoms, those that are furthest from the oxy- 
gen vacancy among these eight titanium atoms, move 
0.09 A toward the strontium vacancy (0.06 A along e 3 , 
and 0.05 A along both e\ and e-i). The final four of these 
eight nearest-neighbor titanium atoms move toward the 
vacancies by 0.11 A (0.05 A along B3, and 0.08 A along 
either e\ or e 2 , depending upon symmetry, with 0.05 A 
along the opposite vector, chosen such that more of each 
displacement in the e\e.2 plane is toward the oxygen va- 
cancy). All other atoms move less than 0.17 A. 



We now investigate correlations between these above 
reconstructionally averaged local displacements and the 
far-field defect-strain tensor. We find surprising results 
in this case of the strontium-oxygen divacancy. While 
the defect-strain tensor and local displacement pattern 
both show expansion on the [001] axis outward from the 
defect, we see a disagreement in sign in the (001) plane 
between the far-field contraction and the expansion of the 
nearest neighbors. This serves as a cautionary note that 
far-field and near-field strain patterns need not be simply 
related and reinforces the importance of providing both 
sets of information for further experimental analysis of 
x-ray scattering signatures. 
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FIGURE 10: Diagonal components of elastic dipole tensor for 
titanium-oxygen divacancy in strontium titanate. Data show 
linear convergence with inverse linear dimension of supercell 
size. 



E. Titanium-oxygen divacancy 



As described above in Section [TTTJ there are six distinct 
orientations for the titanium-oxygen divacancy as defined 
by the V-ri-Vo direction. Moreover, because of the recon- 
struction, there are in fact two distinct classes of sites 
within each possible orientation, as distinguished by the 
rotation state of the octahedron in which the titanium 
sits. Below, we report results for a [100] defect with the 
titanium sitting in an octahedron of positive rotation. 

The elastic dipole tensor for this defect is 



G 



TiO 



14.74 -0.08 
-0.08 18.13 
-4.76 -3.87 




where Figure 10 shows the convergence of the diagonal 
elements of the above tensor. The reconstructional aver- 
age gives 



G 



TiO 



14.74 0.00 0.00 
0.00 19.84 0.00 
0.00 0.00 19.84 



eV, 



with corresponding defect-strain tensor, 



k-TiO 



3.00 
0.00 
0.00 



0.00 
21.68 
0.00 



0.00 
0.00 
21.68 



A 3 



The defect-strain tensor expresses the tendency of this 
defect to expand the crystal in all directions, but primar- 
ily along the directions orthogonal to the V-xi-Vo axis. 
When the above tensor is averaged over all six defect 
orientations (oxygen sites in octahedra surrounding the 
central titanium site, recalling that the Vri-Vo axis is di- 
rection dependent) , the result is a constant tensor with a 
uniform chemical strain per unit concentration of defect 
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FIGURE 11: Local strain pattern for titanium-oxygen diva- 
cancy in the eie2 plane of strontium titanate: titanium atoms 
(♦), oxygen atoms (#), titanium vacancy ( ), oxygen va- 
cancy (O). Atomic displacements exaggerated by a factor of 
three for clarity (->•), displayed for significant in-plane dis- 
placements (> 0.1 A) only. 



of +15.45 A 3 . Such expansion corresponds to a ratio of 
chemical strain e c to the deviation from stoichiometry S 
in SrTii_,503_5 of e c /S — +0.260, indicating a significant 
tendency for the crystal to expand due to the presence 
of titanium-oxygen divacancy. 

Now that our reconstructional averaging has sym- 
metrized the displacement patterns of neighboring atoms, 
it is instructive to examine this set of near-field atomic 

The situation with this diva- 



strains, shown in Figure 11 



cancy is again more complicated than that of earlier iso- 
lated atomic vacancies. There are six atoms in the first 
shell around the divacancy. The five remaining (the sixth 
is itself missing) oxygen atoms that are nearest neighbors 
of the titanium vacancy all move away from that vacancy: 
the one that is furthest from the oxygen vacancy moves 
0.20 A directly away from the vacancy, while the other 
four (which are in the 6263 plane) move 0.23 A away from 
the vacancy (0.13 A along e± toward the oxygen vacancy, 
with the remaining projection of 0.19 A away from the 
divacancy either along 62 or £3 as dictated by the symme- 
try). The sixth atom in this first shell is the sole titanium 
atom that is the nearest neighbor to the oxygen vacancy; 
it moves directly away from the vacancy by 0.13 A. In the 
second shell, we identify four of the oxygen atoms that 
are nearest neighbor to the oxygen vacancy (we already 
counted the other four nearest-neighbor oxygens above, 
"assigning" them to the titanium vacancy) ; these oxygens 
move 0.33 A toward the oxygen vacancy (0.31 A along ei, 



11 



Defect 


tc/5 


Vo 


0.001 


V Sr 


0.030 


V Ti 


0.402 


VsrO 


-0.008 


Vno 


0.260 



TABLE III: Individual ratios of chemical strain e c to stoichi- 
ometric defect deviation S for different defects as calculated 
in Section HVl 

with remaining projection of 0.10 A either along 62 or 83 
as dictated by symmetry). Also, in this second shell, we 
can loosely consider the eight nearest-neighbor strontium 
atoms to the titanium vacancy, all of which move toward 
the titanium vacancy: four of these strontium atoms are 
in a plane that includes the oxygen vacancy, and these 
move 0.22 A toward the titanium vacancy (0.17 A along 
ei, and 0.09 A along both e~2 and I3); the other four 
strontium atoms are on the opposite side of the titanium 
vacancy, and these move significantly toward the vacancy 
by 0.58 A (0.38 A along ei, and 0.31 A along both e2 and 
63). All other atoms move less than 0.20 A. 

We again examine connections between the reconstruc- 
tionally averaged local displacements and the far-field de- 
fect-strain tensor. In this case of the titanium-oxygen di- 
vacancy, the nearest-neighbor atoms to both vacancies 
move outward, showing the same behavior as the de- 
fect-strain tensor. These atoms also demonstrate larger 
movements in those directions orthogonal to the Vri-Vo 
axis, which conforms with the far-field tensor. The next- 
nearest neighbors to both vacancies move inward with 
significant displacements, with those closest to the tita- 
nium vacancy moving by almost three times the amount 
that the nearest-neighbor atoms move. So here we note 
once more that the far-field strain does not appear to 
correlate with the largest magnitude displacement, but 
instead with that of the nearest-neighbor atoms. 



V. DISCUSSION 

The above results for the local strain patterns and de- 
tailed defect tensors are now available for direct compar- 
ison with diffuse x-ray measurements; however, we are 
not aware of any such x-ray data to date. Nonetheless, 
the fully averaged (over both reconstructions and orien- 
tations) defect-strain tensors relate directly to measure- 
ments which are commonly done of chemical strain as a 
function of defect concentration. Table ITnl summaries the 
results from Section |TV| for the ratios of chemical strain e c 
to stoichiometric defect deviation S for all of the defects 
considered in this study. 

Oxygen-vacancy concentration is widely thought to 
serve a crucial role in the properties ofperovskitesj^HII] 
is readily varied but difficult to controlj^ESl and is exper- 
imentally observed to affect chemical strain.^ Moreover, 



Reference state Ae c /S 

Bulk 0.001 

V Sr -0.038 

V Ti -0.142 



TABLE IV: Net ratios of chemical strains to stoichiometric 
defect deviation for oxygen vacancies, referenced against bulk 
and isolated cation vacancies as appropriate. 

cation stoichiometry is also difficult to controP^ISI] and 
so it is uncertain whether, as oxygen vacancies are intro- 
duced into the crystal, such vacancies bind to cation va- 
cancies or form in isolation. In the former case, where the 
oxygen vacancies eventually bind to preexisting cation 
vacancies, the reference configuration should be the crys- 
tal containing the cation vacancy. Hence, it is the dif- 
ference between the chemical strain of the oxygen-cation 
divacancy and that of the isolated cation vacancy that 
describes the change in the crystal lattice as a function 
of varying oxygen-vacancy concentration. In the latter 
case of isolated vacancies, the bulk crystal is in fact the 
system into which these vacancies are introduced, and 
thus the chemical strain as a function of oxygen- vacancy 
concentration is described precisely by that of the iso- 
lated oxygen vacancy in our study. Table |IV| summarizes 
the resulting net chemical strain Ae c versus oxygen-va- 
cancy concentration S. 

For the oxygen vacancy, we have the intriguing result 
that the elastic dipole tensor and corresponding defect- 
strain tensor almost vanish under oricntational averag- 
ing. Thus, very little net effect on the lattice can be 
expected from the presence of isolated oxygen vacancies. 
Moreover, the large anisotropy of the dipole tensor of the 
oxygen vacancy and the ease of introduction and high 
mobility of such vacancies should allow for the control 
of the population and orientation of oxygen vacancies by 
applying external stress. (For instance, at 1% strain, the 
orientational energy differences from Go are 66 meV, or 
about 2.6 times room temperature.) Also, such vacan- 
cies can be expected to tend to shield internal crystalline 
stresses that result from materials processing, a fact po- 
tentially related to the observed difficulties in control- 
ling the oxygen-vacancy concentration during crystalline 
growth. 

One of the earliest sets of available experimental data 
on chemical strain due to oxygen vacancies in strontium 
titanate comes from Yamada and Miller who unfortu- 
nately found a null result. Nonetheless, that null result 
places bounds which, in conjunction with our results, al- 
low some conclusions to be drawn. Yamada and Miller 
varied the oxygen-vacancy concentration over a range 
from nearly zero up to 3.24 x 10 19 cm~ 3 (S = 0.0019 
in SrTiOs-^), stating that "no volume change upon re- 
duction was assumed," due to the experimental uncer- 
tainty of the lattice constant (Aa = 5 x 10~ 4 A) in their 
x-ray diffraction measurements. Their detection limit 
for chemical strain per stoichiometric defect deviation is 



12 



therefore |e c /<S| = \(Aa/a)/5\ < 0.066, where a is the cu- 
bic lattice constant. From Table II VI it is evident that this 



bound is consistent with either isolated oxygen vacancies 
or strontium-oxygen divacancies, but is inconsistent with 
titanium-oxygen divacancies. 

The literature also presents studies of chemical strain 
due to cation non-stoichiometry. Ohnishi ct al. 59 present 
experimental results on the ratio of chemical strain to 
deviation from cation stoichiomctry for samples grown 
by pulsed laser deposition. Specifically, they measure 
the lattice changes for both strontium-rich and stron- 
tium-poor strontium titanate. They associate the stron- 
tium-rich phase with creation of Ruddlesden-Poppei^^ 
planar faults (extra SrO layers) and the strontium-poor 
regime with the presence of strontium vacancies, possi- 
bly bound into strontium-oxygen divacancies. A least- 
squares fit to the results of Ohnishi and coworkers pro- 
vides a value for e c /8 between +0.13 (non-weighted) and 
+0.14 (weighted by reported experimental uncertainty) 
for the strontium-rich phase (Sri+^TiOa) and between 
+0.5 (non-weighted) and +0.8 (weighted by reported 
experimental uncertainty) for the strontium-poor phase 
(Sr!_ 5 Ti0 3 ). 

Our calculations of titanium-vacancy chemical strains 
are not directly relevant to the strontium-rich phase be- 
cause they do not account for Ruddlesden-Popper pla- 
nar faults. On the other hand, our results for strontium 
vacancies, under the interpretation of Ohnishi et al.j^ 
should be d irec tly relevant to their strontium-poor sam- 
ples. Table |m| gives eJS = +0.030 and e c /S = -0.008, 
respectively, for isolated strontium vacancies and bound 
strontium-oxygen divacancies. However, both of these 
are an order of magnitude smaller than the observed 
chemical strain. (We note parenthetically that, while our 
calculations reflect chemical strains for isolated defects 
and Ohnishi and coworkers measure strains for relatively 
high defect densities, Figures [3] and [8] show that trends in 
our data with increasing defect concentration only tend 
to reinforce our conclusions.) We remark that our re- 
sults are consistent with the observation that, within 
a given structural class, the lattice constants of the ti- 
tanates are largely insensitive to the nature of the A-site 
cations. (See data on A 2+ B 4+ 03 perovskites compiled by 
GalassoP^l) From both the results of our calculations and 
this general observation, it seems implausible that simple 
A-site vacancies should produce the measured magnitude 
of chemical strain. Chemical strains of the magnitude 
measured by Ohnishi and coworkers more plausibly arise 
from B-site vacancies or defect complexes associated with 
such vacancies. 

We note that, intriguingly, there is an approximate 
coincidence between the order of magnitude of our calcu- 
lated ratio of chemical strain to stoichiometric deviation 
due to titanium vacancies, |e c /<5| = 0.4, and the observed 
values (+0.5 and +0.8, depending upon weighting of fit) 
for the strontium-poor samples. Given the magnitude of 
the observed lattice expansion and the fact that only tita- 
nium vacancies appear capable of producing an effect of 



this size, we are led to the intriguing conjecture that, per- 
haps, the strontium-poor samples exhibit defects that in- 
clude titanium vacancies and thus a structure more com- 
plex than initially thought (simple strontium vacancies or 
strontium-oxygen divacancies). Clearly, more investiga- 
tion is needed on this point, specifically as to the nature 
of the point defects in the strontium-poor samples grown 
by pulsed laser deposition. 

The above measurements by Ohnishi and coworkers^ 
were performed on strontium titanate samples de- 
posited by the highly energetic process of pulsed laser 
deposition.^ We have learned of recent results by Brooks 
et al.J^ that repeat the above measurements on stron- 
tium titanate samples grown by molecular -beam epitaxy, 
which is a lower-energy deposition process 65 66 and thus 
less prone to the creation of point defects.^ Performing a 
weighted least-squares fit to the new data of Brooks and 
coworkers, we find a value for e c /8 of +0.032+0.019 in the 
dilute limit of the strontium-poor regime. This chemical 
strain ratio, published after our initial calculated predic- 
tion was submitted to this journal, shows good agreement 
with our value of +0.030, and thereby lends convincing 
support for our methodology as well as the applicability 
of empirical shell potentials to calculate reasonable es- 
timates of the chemical strain per stoichiometric defect 
deviation of vacancies. 



VI. SUMMARY AND CONCLUSION 

We have calculated both near- and far-field strains 
for five defects in reconstructed strontium titanate: iso- 
lated oxygen, strontium, and titanium vacancies, as well 
as strontium-oxygen and titanium-oxygen divacancies. 
Given the propensity of the crystal for local fluctuations 
in the reconstruction, we report results both for a partic- 
ular reconstructed state and as averaged over all possible 
local reconstructions. The reconstructionally averaged 
near-field strain results are presented and interpreted in 
terms of the movement of neighboring shells of atoms 
at increasing distances from the vacancy or divacancy. 
We report far-field strain results in terms of both clastic 
dipole tensors, with and without reconstructional aver- 
aging, and associated defect-strain tensors, with recon- 
structional averaging. Anticipating that far-field effects 
will necessarily involve contributions from an ensemble 
of defects, we also present results averaged over all pos- 
sible orientations of the defect within the bulk crystal. 
From these averaged tensors, we extract the resultant ra- 
tio of chemical strain to stoichiometric defect deviation. 
Finally, the combination of local and long-range results 
presented herein will enable determination of x-ray scat- 
tering signatures for comparison with experimental re- 
sults and should further motivate future work on defect 
mechanics, including the influence of externally imposed 
strain (such as in heterostructures) on vacancy popula- 
tions. 

For the oxygen vacancy, we find a highly anisotropic 
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elastic dipolc tensor, with almost perfect cancellation un- 
der orientational averaging. This may be correlated with 
observed difficulties in controlling oxygen concentration 
and lack of observation of effects of oxygen- vacancy con- 
centration on lattice strain. The anisotropy of this tensor 
also suggests that oxygen vacancies may provide a mech- 
anism to shield local internal strains and that applica- 
tion of external stress may allow for control of oxygen 
stoichiometry. From comparison to lattice-strain stud- 
ies, we identify both isolated oxygen vacancies and bound 
strontium-oxygen divacancies as consistent with the ex- 
perimentally observed chemical strain as a function of 
oxygen-vacancy concentration in strontium titanate. 

For cation non-stoichiometry, we find strong indica- 
tions that the point defects in strontium-poor strontium 
titanate samples grown by pulsed laser deposition arc 
not simple strontium vacancies or strontium-oxygen di- 
vacancies, but likely more complicated defect complexes. 
Further, we identify indications that titanium vacancies 



may play a role in these defect complexes. Finally, dur- 
ing the review process, we learned of recent experimental 
data, from strontium titanate films deposited via molec- 
ular-beam epitaxy, that show good agreement with our 
calculated value of the chemical strain associated with 
strontium vacancies. 
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